www.gusucode.com > 基于粒子滤波的故障检测,使用似然函数作为检测函数 > 基于粒子滤波的故障检测,使用似然函数作为检测函数/code/FDI based on SIR likelihood/drawpar.m

    function xnext = drawpar(xprev, stdw, k)
M = length(xprev);
b0=25;
for i=1:M
    %if i<301
     %   b=b0;
    %else
     %   b=b0*0.5;
    %end
xstate = xprev(i)/2 + b0*xprev(i)/(1+xprev(i)^2) +8*cos(1.2*(k-1));
xnext(i,1) = stdw*randn+xstate;
end